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Abstract 

We have performed copper NQR experiments in high temperature supercon- 
ductors YBa2Cu4C>8, YBa2Cu3C>7, and L^-ajSr^CuCU (x=0.12 and £=0.15), 
using the stimulated echo technique which utilizes the rf-pulse sequence 
7r/2 — (t) — 7r/2 — (T — t) — tt/2. The r and T dependences of the stim- 
ulated echo intensity is analyzed by a model that includes the spin-lattice 
relaxation process (Ti-process) and the fluctuating local field due to nuclear 
spin-spin coupling. The model gives quantitative account of the experimental 
results in YBa2Cu40s and YBa2Cu307, using the known values of 1/T\ and 
I/T2G, the gaussian decay rate of the spin echo intensity. The same model 
applied to La2_ x Sr x Cu04 enables us to extract the value of T^g- Our re- 
sults indicate that T1T/T2G is independent of temperature, implying that the 
dynamic exponent is one in L^-zSr^CuC^. 
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I. INTRODUCTION 



Anomalous spin dynamics in high temperature superconductors continues to be a contro- 
versial issueHH. In all high-T c materials, the dynamic susceptibility x(Qj u ) a t low frequencies 
is peaked near or at the antiferromagnetic wave vector Q = (tt, it) above T c . The nuclear 
spin-lattice relaxation rate (1/Ti) at the planar 63 Cu sites measured by NMR/NQR has 
been used extensively to probe the low frequency spin fluctuations through the following 
relation! 



1 = T^B 

TiT & 



E^te)V(g,"»)M» (i.i) 



where 7 n is the nuclear gyromagnetic ratio, u n is the nuclear resonance frequency, and F±(q) 
(F||) is the Fourier component of the hyperfme coupling perpendicular (parallel) to the crys- 
talline c-axis, which is assumed to be the quantization axis of nuclear spins as is the case 
for NQR. It was found that 1/{T]T) shows contrasting temperature dependence for different 
regions of doped hole concentration. For sufficiently high hole concentration, 1/(T\T) keeps 
increasing as temperature is lowered down to T c . In contrast, in many underdoped materials, 
where the hole concentration is smaller than the optimum value for highest T c , 1/(T]T) shows 
a broad maximum at a temperature T s , which is much higher than T c , and decrease steeply 
at lower temperatures. This suggests that a pseudo gap opens in the spin excitation spec- 
trum near (71*, 7r) above tJI. As the hole concentration increases, T s decreases. The pseudo 
spin-gap phenomenon is observed in various underdoped materials, such as YBa 2 Cu 3 66 3E, 
YBa 2 Cu 4 8 l, LaBa 2 Cu 3 7 J and HgBa 2 Ca n _ 1 Cu n 2n+2+(5 (ri = 1,2)11. 

However, one typical high-T c system, La 2 _ x Sr x Cu04 (LSCO) does not show clear pseudo 
spin-gap behaviortlJ. Instead, recent precise measurements of l/(T\T) up to sooidlHi 
revealed a modest change of behavior at a temperature T*. Above T*, 1/(T\T) follows 
a Curie- Weiss temperature dependence, 1/(T + 0). Below T* , \/(T\T) keeps increasing 
but is slightly suppressed from the extrapolation of the high temperature Curie- Weiss law. 
Dependence of T* on the hole concentration is similar to that of T s in other underdoped 
high-T c materials. At present, origin of this crossover temperature or the reason for the 
different behavior between LSCO and other high-T c materials is not well understood. 

NQR experiments in LSCO are complicated due to large inhomogeneous line width of 
a few MHz. Since Cu spins are excited over a few hundred kHz by typical rf-pulses, only 
a small portion of the entire Cu spins contributes to the NQR signal at one frequency. In 
our previous work, it was found that 1/Ti measured by Cu-NQR is not uniform over the 
broad spectrum. The frequency dependence of 1/Ti is more pronounced in more underdoped 
materialsEl, indicating substantial inhomogeneity of the microscopic electronic state. One 
concern is that the measured value of 1/Ti might be influenced by the spectral diffusion 
process, in which Zeeman energy Ti^ n I z H of the excited nuclear spins diffuses in frequency 
over inhomogeneously broadened spectrum by the transverse component of nuclear spin-spin 
coupling of the form al i+ lj_. If this is the case, the measured value of 1/Ti would be larger 
than the true relaxation rate due to electronic spin fluctuations. 

Another issue is that one cannot obtain the strength of nuclear spin-spin coupling from 
measurements of the spin echo decay rate for such a broad spectrum. In cuprates, Cu 
nuclear spins are coupled dominantly by the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
type indirect interaction mediated by electronic spin excitations, 
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H = hJ2 {aptzljz + o?/2 (Ii+Ij- + h-I j+ )} , 



;i.2) 



where the coupling constant is given in terms of the g-dependent static susceptibility as 

1 



-, l J 



H9V>b)' 



E \ F M)\ 2 x{q) exp(-iqRij). 



1.3) 



The large anisotropy F\\(q)/F±(q) ~ 6 for q ~ (tt,tt) makes an much larger than a±. Then 
the Gaussian time constant of spin echo decay 1/T 2 g in 63 Cu NQR experiments is given 



T, 



2G 



0.69 
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where 0.69 is the natural abundance of 63 Cu nucleus. Measurements of 1/T 2 q have pro- 
vided valuable information regarding the enhancement of the static susceptibility near the 
antiferromagnetic wave vector in many high-T c materials such as YBa2Cu30"7_AiY1237, 
T c ~ 90^)00, YBa 2 Cu 3 6 . 62 (T c ~ 60^)0, and YBa 2 Cu 4 8 (Y1248, T c ~ 80AT)S How- 
ever, Eq. (|1.4|) assumes that the nuclear spins over the entire spectrum are flipped by the 
7r-pulse, a condition that is not satisfied in LSCO. 

In this paper, we report the results of stimulated ec hc§ decay measurements in 
YBa 2 Cu40s, YBa 2 Cu3C>7, and La 2 _ x Sr :c Cu04 (x—0.12 and x=0.15), using NQR at zero 
magnetic field. Stimulated echo sequence utilizes three tt/2 pulses separated in time by r 
and T — t. The decay of stimulated echo intensity as a function of T can be caused not 
only by the spin-lattice relaxation process but also by the spectral diffusionS as well as the 
nuclear spin-spin coupling0. The dependence on r and T of the stimulated echo intensity 
is compared with the calculation based on a model which takes account of the spin-lattice 
relaxation process and the nuclear spin-spin coupling along the c-axis. The effects of spin- 
spin coupling perpendicular to the c-axis, which is much smaller than the coupling along the 
c-axis but responsible for the spectral diffusion, is neglected as we discuss in detail below. 
For YBa 2 Cu40§ and YBa 2 Cu30"7, where the values of 1/Xi and 1/T 2 g are known, the cal- 
culation shows good quantitative agreement with the experimental data. The same model 
applied to LSCO enables us to extract the value of 1/T 2 g- Our results indicate that T\T/T 2 g 
is independent of temperature, implying that the dynamical exponent is one in LSCO. 



II. SAMPLE PREPARATIONS 



We used polycrystalline materials synthesized by the usual solid state reaction. For 
Y1248 material, we used hot-isostatic-pressing (HIP) technique. All samples were confirmed 
to be single phase by powder X-ray diffraction measurements. The superconducting tran- 
sition temperatures (T c ) are 90K (Y1237), 80K (Y1248), 38K (LSCO x = 0.15) and 32K 
(LSCO x = 0.12), respectively. 
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III. STIMULATED ECHO METHOD 



We briefly discuss the stimulated echo method in comparison with the ordinary spin echo 
and inversion recovery methods. In ordinary spin echo experiments, two rf-pulses (tt/2 and 
7r pulses) are applied with time separation r. The r dependence of spin echo intensity is 
described by two parameters, 1/T 2 l and 1/T 2 g as0, 

M(2r) = M exp(-|^)exp( — . (3.1) 



2L 



The Lorenzian component 1/T 2 l represents the contribution from the spin-lattice relaxation 
process^. The Gaussian component 1/T 2 q is due to nuclear spin-spin coupling and given 
by Eq. (0). 

The pulse sequence for the stimulated echo experiments are described as (vr/2 — r — n/2 — 
(T — r) — vr/2), which is shown in Fig. |l](b) and compared with the conventional inversion 
recovery sequence (tt — T — n/2 — t — ir) in Fig. |T](a) . The figures below the time chart of 
rf-pulses show the corresponding state of nuclear spins in the rotating frame. In the inversion 
recovery sequence (a), the nuclear magnetization along the ^-direction initially at the thermal 
equilibrium is inverted by the first 7r pulse. The ^-component of the magnetization M z , which 
then begins to recover toward the thermal equilibrium, is measured by the intensity of the 
spin-echo signal formed by the second and the third pulses. In case of NQR of 63 Cu nuclei 
(/ = 3/2), the recovery of spin echo intensity is described by an exponential function if the 
effect of spectral diffusion is neglected, 



M(T) = M 



1 — 2 exp 



3T 



(3.2) 



where T\ is the time constant for high field NMR experiments. 

In the stimulated echo sequence (b), the initial magnetization is flipped into the y- 
direction by the first 7r/2 pulse. Then it dephases in the xj/-plane by the inhomogeneous 
distribution of the resonance frequency, as well as by the nuclear spin-spin coupling. The 
magnetization spread in the x|/-plane is flipped into the xz-plane by the second 7r/2 pulse. 
The value of M z immediately after the second pulse is given by the ensemble average of 
— cos0(r), where 0(r) = 7 n Jq h{t)dt is the accumulated phase in the xy-plane of individual 
spins subject to the local field h(t) during the period between the first and second pulses. 
If the local field distribution is dominated by the static inhomogeneity of the resonance 
frequency, M z should be oscillating as — cos(ru;) as a function of frequency over the spectrum 
as shown in Fig. |2|. If M z does not change between the second and the third pulse, the 
oscillating magnetization distribution along the z-direction is flipped into the ^/-direction by 
the third pulse, which then refocuses along the — y-direction at T + r, forming the stimulated 
echo. The echo intensity is given by the ensemble average of (1/2) cos0(T + r), where 

<P(T + r) = 7n f h(t)dt - ln j T+T h(t)dt (3.3) 



JT 



is the accumulated phase at T + r. The accumulated phase is zero if h(t) is static. 

Temporal variation of M z or h(t) causes the intensity of stimulated echo to decay as 
a function of T. Change of M z can be caused either by the spin-lattice relaxation process 
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(Ti-process) or by the spectral diffusion. By the Ti-process M z recovers towards the uniform 
thermal equilibrium magnetization. The stimulated echo decay by the ^-process is simply 
given by exp(— 3T/Ti) similar to the inversion recovery. 

The spectral diffusion can be caused by the transverse component of the nuclear spin- 
spin coupling, the second term of Eq. ( |1.2| ), which describes the mutual flip of two nuclear 
spins. We discuss this in more detail in the next section. The oscillating distribution of M z 
will be averaged out by the spectral diffusion, leading to reduction of the stimulated echo 
intensity. 

Stimulated echo decay is also caused by the temporal fluctuation of h(t), since the two 
terms in the accumulated phase Eq. ( |3.3| ) then do not cancel exactly. In our NQR exper- 
iments, distribution of the local field is largely due to inhomogeneity of the electric field 
gradient, which is static. However, the longitudinal component of the spin-spin coupling, 
the first term of Eq. (|1.2| ), gives rise to an additional source of the local field. This local 
field at i-th site is given as 

M*)=-£°y j i*(*)> (3-4) 

which fluctuates via the Ti-process. 



IV. EXPERIMENTAL RESULTS AND ANALYSIS 

The stimulated echo intensity M was measured as a function of T for several different 
values of r. The experimental results are shown in Figs. (| for Y-based compounds and 
Fig. for LSCO x = 0.12 sample. Generally, the echo decay curve M(T) deviates from a 
single exponential function exp(— 3T/Ti). The deviation is not clearly noticed for small r but 
becomes significant for larger r. The decay rate d\aM{T)/dT shows strong r dependence 
for small T, where it is larger for larger r. In contrast, the decay rate depends less on r for 
larger T region. How strongly the decay curve changes with r depends on temperature as 
well as materials. As shown in Figs. |3], f|, and the r dependence of stimulated echo decay 
curves become negligible at higher temperatures. At lower temperatures, the r dependence 
of decay curves is most significant in Y1248 sample. 

As mentioned in the previous section, we have to consider three distinct processes for the 
stimulated echo decay, namely, the Ti-process, spectral diffusion, and local field fluctuations. 
Of these, contribution from the Tx-process is easily separated by dividing the experimental 
data of stimulated echo decay by exp(— 3T/Ti). 

Portis has analyzed the spectral diffusion process in electron paramagnetic resonance of 
F-centers, where diffusion occurs via the dipolar interaction between electron spins within 
the spectra broadened by the hyperfine interaction with surrounding nuclear spins^. We 
apply this analysis to the present problem. First, the inverse life time of the z-component 
of a nuclear spin due to mutual spin flip of Ii and a nearby spin Ij described by the second 
term of Eq. ( |1.2|) is given byH 

i = 0.69x^/M£(aY) 2 , (4.1) 

s j 
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where f{oj) is the NQR spectral shape function normalized so that / f(u)du = 1 and we 
assume that the resonance frequencies of two neighboring nuclear spins are uncorrelated. 
Although the total Zeeman energy has to be conserved, the uncertainty principle allows the 
resonance frequencies of Ii and Ij to be different by an amount comparable to l/r s . Thus 
the mutual spin flip process describes a one dimensional random walk of magnetic particles 
along the frequency axis, where both the hopping rate and the hopping distance are given by 
1/t s . Therefore, the time evolution of the macroscopic magnetization distribution is given 
by a diffusion equation with the diffusion constant D w = We can estimate l/r s as 

follows. From Eq. (|1.4|), 



r s \J-2gJ Y*tj(a(r 

where E^^af) 2 /^^) 2 ~ {F±{Q) / F\\(Q)V ~ 1/40 in high-T c cuprates. Among the 
materials we studied, r s should be shortest in Y1248 because of small NQR line width of 
200kHz (f(u) = (4vr x 10 5 )" 1 ) and relatively large 1/T 2G ~ 3 x 10 4 sec" 1 near T c . From 
these numbers, l/r s is estimated to be 56 sec -1 . 

In principle, such spectral diffusion could average out the oscillation of the magnetization 
distribution along the frequency axis with the period 27r/r in the stimulated echo experi- 
ments (Fig. |]). The characteristic time for this to occur is given by (tt/t) 2 / D^- Even for 
the longest value of r = 5 x 10 _5 sec in our experiments, this characteristic time is of the 
order of 10 4 sec, which is many orders of magnitude larger than the range of T in our exper- 
iments (less than 10 _2 sec). Therefore, we conclude that the effect of the spectral diffusion 
is negligible. 

We now consider the effect of fluctuations of the local field. Recchia et al. have considered 
the same problem for the ordinary spin echo decayl 4 !. They assumed a Gaussian distribution 
for the accumulated phase and derived an analytic expression for the spin echo decay in terms 
of the correlation function of the local field. Curro et al. applied this results to the case 
where the local field is given by Eq. (|3.4| ) and I z of a neighboring nuclear spin is fluctuating 
via the T\ processi. This approach can be applied to the stimulated echo decay with only 
minor modification. If we assume a Gaussian distribution for the accumulated phase defined 
by Eq. (|3.3|) , the stimulated echo intensity is given by 

(ico^(T + r)) = iexp(-MI±l)!>j. (4.3) 

As discussed in the appendix, (<p(T + r) 2 ) is the ensemble average of the contribution from 
individual neighboring nuclear spins. Each of them has to be distinguished whether it is a 
like nucleus, which is on resonance, or an unlike nucleus, which is off resonance. Thus if 
there are n types of unlike nuclei, we can write 

n 

(0(T + r) 2 ) = P ((puUT + rf) + £ Pi(4> unlike ,i(T + r) 2 ), (4.4) 

i=i 

where Pq and Pi are the abundance of the like nuclei and z-th unlike nuclei, respectively. In 
the appendix, (4>uke(T '+ t) 2 ) and ((f> U nUke,i(T + t) 2 ) are calculated as functions of T 2 g of 63 Cu 
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nuclei as defined in Eq. (|1.4j) and Ti's of like and unlike nuclei. In Fig. [g, T dependence of 
(4>Uke(T + t) 2 ) and {4> U niike(T + r) 2 ) are plotted for different values of T 2 g/Ti and r/Ti. 

To summarize our analysis, we expect that the T and r dependence of the stimulated 
echo intensity M(T, r) is given by the product of the two decay process, i.e. the Tx-process 
and local field fluctuations, 

MC r, T ) = «b«p(-f)«p(-^±^>). (4.5) 



V. DISCUSSIONS 

In this section we compare the experimental results with the calculation of Eq. ( ,5| ). 
In order to examine the validity of our approach to analyze the local field fluctuations, we 
separate the trivial Tx-process. In the following, we compare the stimulated echo intensity 
divided by exp(— 3T/Ti), which we call the corrected intensity M corr (T,r), with the calcu- 
lated results of exp(— (<p(T + r) 2 )/2). The results of comparison are shown in Figs. [7|, ||, 
and |9]. 



A. Y-based compounds 



We first discuss the results on underdoped YE^Ou/jOg and optimally doped 
YBa 2 Cu307_,5 (T C =90K). The full width at half maximum of 63 Cu-NQR spectra is about 
200 kHz in these samples. Since almost all 63 Cu nuclear spins are excited by rf-pulses for 
such a narrow spectrum, all 63 Cu spins are like nuclei and 65 Cu spins are the only unlike 
nuclei. Thus we set P = 0.69 and Pi = 0.31 in Eq. (fyj). In Fig. |7|, the data for M corr and 
the calculated stimulated echo decay curves are compared in Fig. ^ for Y1248 and Y1237 
compounds. For calculation, we used the value of 7\ obtained from NQR inversion recov- 
ery measurements and the data of T 2G obtained by spin echo decay measurement M. The 
overall magnitude M is the only adjustable parameter in this comparison. The agreement 
between the experimental data and calculation is quite good. The calculation reproduces 
the observed change of behavior for different values of r, in particular, the non-exponential 
decay for large values of r is well explained by our model. Thus we believe that our model 
is valid for these materials. 



B. La2_ x Sr z Cu04 



As we mentioned before, the Cu-NQR spectra of LSCO show significant inhomogeneous 
broadening, whose full width at half maximum (FWHM) is about 2MHz. Thus it is not 
possible to obtain 1/T 2 g as defined in Eq. (|1.4p. On the other hand, Walstedt et al. measured 
the spin echo decay rate using the I g = 1/2 -1/2 center line of the 63 Cu NMR in the 
oriented powder sample in high magnetic field!!. The FWHM of the NMR center line in 
LSCO is typically about 500Oe. Although this is about 4 times narrower than that of Cu 
NQR, the whole region of the spectra can not be excited by rf-pulses. Moreover, the NMR 
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spectra includes both the broad background from unoriented portion of the sample and the 
resonance from the so called .B-site, which may affect the spin echo decay measurements. 
Thus it would be important to obtain l/T 2G by a different method. 

Here we assume that our model for the stimulated echo decay is valid also for LSCO and 
extract the value of I/I2G by fitting the data of M corr to exp(— (<p(T + r) 2 )/2) with T 2G as 
a fitting parameter. The experimental data and fitted curves for LSCO are shown in Fig. || 
We used the value of T\ determined from the inversion recovery technique. 

There is minor ambiguity in this procedure. Since the 63 Cu NQR spectrum is so broad, 
only a fraction of 63 Cu spins is considered to be like nuclei. However, this is not a serious 
problem, since (4>uke{T + r) 2 ) and (4> U nUke(T + r) 2 ) are nearly identical for the range of 
parameters of our interest if the values of 1/Ti and 7™ are the same, as shown in Figs. |6] and 
|9]. Based on the ratio between the NQR line width and the magnitude of the rf magnetic 
field of our experiments (~ 200kHz), we chose rather arbitrarily that 10% of 63 Cu are like 
nuclei and 90% are unlike nuclei. The obtained value of 1/T 2 g hardly depends on these 
numbers. We show in Fig. |9| two fits obtained by assuming that 100% of 63 Cu are like or 
unlike nuclei. The fitted values of T 2G differ only by 6%. Our analysis to determine T 2 g is 
limited below 150K for x = 0.15 and below 250K for x = 0.12. At higher temperatures, 
where l/T 2G gets smaller, T dependence of M corr becomes too weak to determine 1/T 2G 
reliably. The values of T 2G thus derived are plotted against temperature for two samples, 
x = 0.12 and x = 0.15 in Fig- El- 1/^2G increase steeply as temperature decreases, similar 
to other high-T c materialsEJ'EMa. The results obtained by Walstedt et al. for x = 0.15 
between T = 100K and 300K using high field NMR are also shown by the dashed line, 
after multiplied by a/2 in order to take account of the difference between NQR and NMR 
measurement s@H 

Since the antiferromagnetic correlation is quite strong in high-T c cuprates, we generally 
expect that the dynamic spin correlations obey some kind of scaling relation, which relates 
temperature dependence of various magnetic quantities to those of the antiferromagnetic 
correlation length £ and the characteristic energy of the antiferromagnetic spin fluctuations 
u s f. These two quantities are related by the dynamical exponent z as uj s f oc £~ z . The data 
of 1/T 2 g and \jT{T have been used to extract the value of z in various high-T c cuprates. It 
has been shown that if xiQ) °^ £ 2 > 

1 1 

?2G 



oc^ s/ ,— oc£. (5.1) 
J- 1 J- i-% 



Therefore, we expect T{I '/Tg G to be temperature independent^. While z — 1 is expected 
at the quantum critical point (the boundary between Neel ordered states and disordered 
states) in two dimensional quantum spin systemil, antiferromagnetic spin fluctuations in 
itinerant electron systems is generally believed to be described by z = 2. 

Experimentally, the underdoped materials of the Y-based system, YBa2Cu3 06.63 and 
YBa2Cu408 shows the z=l behavior above T s . Early data on the optimally doped 
YBa2Cu306.g are consistent with z=2@H. However, recent 17 0-NMR experiments by Keren 
et al. pointed out importance of fluctuations of Cu nuclear spins due to both spin-lattice 
and mutual flip process in the analysis of the spin echo decay datcS They obtained z=l 
after correcting these effectsS. 

In fig.|IT, we plot the ratio T{I '/T^q for £=0.12 and 0.15 against temperature with both 
z=l and z—2. Apparently, the plot for z=l is more temperature independent than the plot 
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for z = 2. Thus our results suggest the quantum critical scaling for the magnetic excitations 
in underdoped LSCO. The results obtained by Walstedt et al. for x = 0.15 between 100K 
and 300K are also consistent with z = lB To summarize our results on LSCO, consistent 
account of the stimulated echo decay data is given by our model using the values of l/Ti 
measured by the inversion recovery technique, which show no clear signature of pseudo spin- 
gap above T c . The quantum critical scaling we found for LSCO is also consistent with the 
lack of spin-gap, since such scaling breaks down in many of the underdoped cuprates below 
T s . Why the spin gap is absent in LSCO is not fully understood. However, we speculate that 
it may be closely related the instability toward the simultaneous spin and charge ordering 
found experimentally in LSCO§. 

VI. CONCLUSION 

We have applied the stimulated echo technique to the high-T c superconductors by 63 Cu 
NQR. Our model assumes that the stimulated echo decay is caused by two dominant pro- 
cesses, the longitudinal T\ relaxations and fluctuation of local field produced by longitudinal 
component of the nuclear spin-spin coupling. We argue that effects of spectral diffusion due 
to transverse component of the spin-spin coupling is negligible. We confirm that our model 
indeed gives good quantitative account of the experimental data for Y-based materials. We 
used the same model to deduce the value of 1/T 2 g from the data of stimulated echo decay 
in LSCO. We found that I/T2G increases monotonously with decreasing temperature. The 
results show that the ratio T{T '/T^q with z=l is approximately independent of temperature, 
suggesting the quantum critical scaling in underdoped LSCO. 

ACKNOWLEDGMENTS 

We are grateful to Prof. H. Yasuoka and Dr. Y. Itoh for suggesting stimulated echo 
experiments on high-T c superconductors. One of the authors (S. F.) thanks to S. Onoda 
for stimulating discussions. This work is supported by the Grant in Aid of the Ministry of 
Education, Sports and Culture. One of the authors (S. F.) is supported by JSPS Research 
Fellowship for Young Scientists. 

APPENDIX: 

We calculate the second moment of the accumulated phase defined in Eq. (|3~3|), taking 
account of the fluctuation of the local field defined in Eq. (|3.4j) due to the T\ process. We 
follow the approach formulated by Curro et a/.i. Since the values of Ij Z is jumping among 
the four eigenvalues ±3/2, ±1/2, the local field is written as 

^W—EE^^W, (ai) 

where Pj m is a probability that Ij z takes the eigenvalue m. From Eqs. ( |A"I| ) and fl3.3|), the 
second moment of the accumulated phase is given as 
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where we have assumed no correlation between different nuclear spins. 
The time evolution of p m (t) is governed by the rate equation, 

dp m (t) „ 

— — = ^ VK m „p n (t) 



(A2) 



(A3) 



W 



2Ti 



^-3 3 ^ 
3-740 
4-73 
V 3 -3 / 



(A4) 



We can write {p m {t)Pm'(t')) = P mm '{\t — £'|)/4, where P mm /(t) is the probability that if a 
nuclear spin is in the state m at the time t under the condition that it was in the state ml 
at t = 0. Pmm'(t) is equal to p m (t) obtained by solving the rate equation, Eq. (|A3|) , with 
the initial condition that p m > = 1 at t = 0. For instance, 



Pz 3 

2 2 



, , 1 9 , t . 1 , 3t, 1 . 6t. 
(*) = 4 + ^ ex P("^) + 4 6XP( -T1 } + 20 eXp( "T\ 



(A5) 



We now have to consider the effects of three n/2 pulses , which redistribute the popula- 
tions of like nuclear spins, which are on resonance. The rate equation, Eq. (|A3|) , should be 
combined with the boundary conditions at t — r, 



P3/2<r + ) =Pi/ 2 (r + ) = ^[p 3 /2(T 



) + 



Pi/a(T )]• 



(A6) 



P-3/2(r + ) =^i/ 2 (r + ) = 2^-3/2^ )+p_i/ 2 (r )]. 



(A7) 



and the similar conditions at t = T + t. The rf-pulses have no effect for unlike nuclear spins, 
which are off resonance. 

Eq. ([AJD shows that (0(T+r) 2 ) is the sum of the contribution from individual neighboring 
nuclear spins. Therefore it is written as Eq. ( f4.4| ) according to the distribution of different 
types of like and unlike nuclei. We define 4>nke(T + r) to be the accumulated phase for the 
hypothetical case where like nuclei occupy all the neighboring sites. Likewise (f) U nUke(T + r ) 
is defined to be the accumulated phase when all neighboring sites are occupied by z-th unlike 
nuclei. They are calculated as, 



2Ti 
0.69 

lexp 



2 r 





2t 


[;{ 




T — 




Zi J 




(A8) 
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and 



unlike. 

i(T + T)) 



/2T U 



2 x 0.69 \ T 2G 
. T + r 

— exp 



+ 2 exp 




T 



it 



2 + 2 exp - 



T 



(A9) 



Here 1/T 2 g is the 63 Cu NQR spin echo decay rate as defined by Eq. (|1.4f ), 1/Tio and 
1/Tu is the spin-lattice relaxation rate of like and 2-th unlike nuclei, and 70 and 7^ is 
the gyromagnetic ratio of like ( 63 Cu) and i-th unlike nuclei respectively. The result for 
j(T + r))) has obtained by Recchia et al. previously^. 
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FIGURES 



FIG. 1. Time charts of three rf-pulses and corresponding states of nuclear spins in the rotating 
frame for (a) Inversion recovery and (b) stimulated echo decay sequences. 

FIG. 2. Distribution of M z in the inhomogeneously broadened spectrum in stimulated echo 
experiments after the second pulse. Typical values of the t w (tt/2 pulse width) and r are 1.5/is 
and 30fis. 

FIG. 3. Stimulated echo intensity is plotted against T for various values of r for YBa2Cu40s 
at 100K. 

FIG. 4. Stimulated echo intensity for YBa 2 Cu 3 7 at 100K. 

FIG. 5. Stimulated echo intensity for La2- x Sr :E Cu04 x = 0.12 sample at 100K (above) and at 
300K (below). 

FIG. 6. The second moment of the accumulated phase (4> 2 ) as a function of T/T\ for different 
values of T2G/T1 and t/T\ for (a)like-spin coupling and (b)unlike-spin coupling. 

FIG. 7. The experimental stimulated echo decay, M corr (dots) is compared with the calculations 
(lines) for Y1248 and Y1237 at 100K. We used the value of parameter, Tig given in ref.18. 

FIG. 8. The experimental stimulated echo decay, M corr (dots) and fitted curves based on Eq. 

4.5. 

FIG. 9. The experimental data of M corr for x = 0.15 at 100K (dots) are compared with two 
fitted curves, assuming 100% like spin coupling (solid curve) or 100% unlike spin coupling (dotted 
curve). The fitted values of T2G are 25.0[/is] for like-spin coupling and 26.7[/xs] for unlike-spin 
coupling. 

FIG. 10. Temperature dependence of \jTiQ for LSCO x = 0.12 and x = 0.15 samples. The 
dotted line shows the results by Walstedt et al. for x = 0.15 obtained from the spin echo decay 
measurements using high field NMR. 

FIG. 11. The plot of T{I '/T* G (z=l or 2) in LSCO x = 0.12 and 0.15 against temperature. 
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